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ABSTRACT 

We derived the absolute proper motion (PM) of the old, solar-metallicity Galactic open cluster M 67 using observations collected 
with CFHT (1997) and with LBT (2007). About 50 galaxies with relatively sharp nuclei allow us to determine the absolute PM of the 
cluster. We find (pi„cos6, yUa)j2ooo.o = (-9.6 ± 1.1, -3.7 ± 0.8) masyr"'. By adopting a line-of-sight velocity of 33.78 ±0.18 kms"', 
and assuming a distance of 81 5±50pc, we explore the influence of the Galactic potential, with and without the bar and/or spiral arms, 
on the galactic orbit of the cluster. 

Key words. Astrometry - Galaxy: open cluster and associations: NGC 2682 (M 67) ~ Galaxy: kinematics and dynamics 
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1. Introduction 

The solar-metallicity Galactic open cluster M 67 (NGC 2682) is 
among the most-studied Galactic open clusters. Still, its absolute 
proper motion (PM) remains poorly constrained. 

We applied for the first time on ground-based multi- 
epoch CCD wide-field images the PM techniques developed in 
Anderson et al. (2006, Paper I), Yadav et al. (2008, Paper II), 
Bellini et al. (2009, Paper III) to define an absolute reference 
frame using background faint galaxies. In the next section we 
describe the data set and the measurements, and a final section 
is dedicated to the study of the orbit of M 67 within the Galaxy 
under different assumptions for the Galactic potential. 



2. Observations, data reduction, proper motions 

Two data sets - collected with two different telescopes and at 
two different epochs - were used to measure the PM of objects 
in the field of M 67. 

As the first epoch (hereafter, epoch 1) we employed images 
taken on Jan. 10-13, 1997, at the CFHT 3.6m telescope. These 
images were first presented by Richer et al. ( I1998I I. We took only 



* Based on data acquired with the Large Binocular Telescope 
(LBT) at Mt. Graham, Arizona, under the Commissioning of the 
Large Binocular Blue Camera. The LBT is an international collab- 
oration among institutions in the United States, Italy and Germany. 
LBT Corporation partners are the University of Arizona on behalf 
of the Arizona university system; Istituto Nazionale di Astrofisica, 
Italy; LBT Beteiligungsgesellschaft, Germany, representing the Max- 
Planck Society, the Astrophysical Institute Potsdam, and Heidelberg 
University; the Ohio State University, and the Research Corporation, 
on behalf of the University of Notre Dame, University of Minnesota 
and University of Virginia; and on observations obtained at the Canada- 
France-Hawaii Telescope (CFHT), which is operated by the National 
Research Council of Canada, the Institut National des Sciences de 
rUnivers of the Centre National de la Recherche Scientifique of France, 
and the University of Hawaii. 

** Visiting PhD Student at STScI under the "2008 graduate research 
assistantship " program. 



a subsample of this data set, specifically 15 exposures of 1200 s 
in the V filter, with a median value for seeing and airmass of 1" 
and 1.2 respectively. Each image was collected with the UH8K 
camera (8 CCDs, 2Kx4K pixels each, with an average scale of 
210 mas pixeP') covering a field-of-view (FoV) of ~29'x29'. 

The second-epoch (epoch 2) data were collected between 
Feb. 16 and Mar. 18 2007 and consist of 56 images of 180 s 
exposures in the B-band filter and 42 exposures of 100-110 s 
in the V-band, obtained with the LBC-blue camera (4 CCDs, 
2Kx4.5K pixels each, FoV of ~24'x26')- A large dither pattern 
(~30% of the FoV) was employed for both filters. Median see- 
ing and airmass were 1" and 1.1 for the V, and r.'3 and 1.1 for 
the B images. We selected only B images for finding objects, 
while those in V were taken to provide all the astrometric in- 
formation. A more extensive description of the LBT data set is 
given in three companion papers (Bellini & Bedin Isubmittedl 
Bellini et al. 120101 and Bellini et al. in prep.). Procedures and 
algorithms used to derive the list of objects, star positions and 
fluxes, and PMs are those explained in great detail in Paper I and 
in Anderson et al. (2008, A08). Below we briefly describe these 
reduction procedures, which were organized in three steps. 

In the first step, we employed the software described in 
Paper I to obtain PSFs, star positions and fluxes in each chip 
of each exposure and for the best sources (bright, isolated, with 
a stellar profile). We corrected LBC raw positions for geometric 
distortion (GD), as described in Bellini & Bedin jsubmittedl l. 
and we used these corrected positions to register all of the 
LBC single-chip images into a common distortion-free refer- 
ence frame (the master frame) by means of linear transforma- 
tions. With this master frame we derived the GD coiTection also 
for the UH8K positions, using the same technique (used also in 
Bellini & Bedin |2009] l. 

The second step consisted of producing the list of objects 
(the target list). We took every single local maximum detected 
within each individual chip of each LBC B image to build a 
peak-map image (A08). The peak map consists of an image with 
the same pixel scale of the master frame, where we added 1 to 
a pixel count each time a local maximum, measured in a given 
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Fig. 1. Vector-point diagram of the absolute PMs, in equatorial 
coordinates (J2000.0). Small crosses are M 67 members (within 
the solid circle), filled circles are the reference galaxies, 47 of 
which (marked with red open circles) were taken to compute the 
centroid of the galaxy distribution. 



image, fell within that pixel (once transformed with the afore- 
mentioned linear transformations). A 3x3 box-car filter is ap- 
plied to the peak map, and in our final list we considered as sig- 
nificant any object 3.2cr (where cr is the rms of the peak-map 
background) above the local surrounding. Our list contains ob- 
jects generating a local maximum in at least 40% of the images 
covering the patch of sky where the maximum is found. We ad- 
ditionally required a minimum coverage of 10 B images. 

With the B images we generated the target list because they 
are more numerous and have a lower background compared to V 
images. In addition, large dithers helped us to produce a list of 
solid detections even for the faintest sources. We further purged 
our list from excessively faint sources, keeping only objects with 
at least 100 DN above the local background and which showed 
up in at least one observation for each of the two V-filter epochs. 
We also excluded objects within chip # 2 and 4 of the UH8K 
camera, because these chips are highly affected by charge trans- 
fer inefficiency. As a consequence, the useful FoV is reduced by 
25%. 

In the last step we derived PMs (see Sect. 7 of Paper I for 
more details): We measured each object in the target list in each 
chip of each V exposure where it could be found, using PSF- 
fitting to get a chip-based flux and a GD-corrected position. Then 
we organized the images in pairs of one image from the first 
and one from the second epoch. For each object in each pair 
we computed the displacement (in the reference frame) between 
where epoch 1 predicts the object position in epoch 2, and the 
actually-observed position in epoch 2. Multiple measurements 
of displacements for the same object were then used to compute 
average displacements and rms. 

It is clear that to make these predictions, we needed a set of 
objects as a reference to compute positional transformations be- 
tween the two epochs for each source. The cluster members of 




4000 5000 6000 7000 8000 9000 10000 



Y [liHC Tiixel] 





m „ - --r 



X [LBC pixel] 

Fig. 2. From top to bottom, fig and jja cos 6 as a function of 
X and Y. When PMs are plotted as a function of their posi- 
tion in the master frame (X,Y), no clear systematic errors ap- 
pear. Note that two chips of UH8K are not used (Y>7000, for 
X<5500 & X>9500), and that the cluster center is at location 
(X,Y)^(7500,7000). These two eff'ects both create gaps in the 
galaxy spatial distribution. 



M 67 were a natural choice, as their internal motion is within our 
measurement errors (~0.2 mas yr ' Girard et al. I1989I I. provid- 
ing an almost rigid reference system with the common systemic 
motion of the cluster. 

We initially identified cluster members according to their lo- 
cation in the V \s. B - V color-magnitude diagram. We took 
only main sequence (MS) stars to transform each exposure into 
the master frame, because it spans a narrow range in coloo By 
predominantly using cluster members, we ensured the PMs to be 
measured relative to the bulk motion of the cluster We iteratively 
removed from the member list those objects with a field-type 
motion (i.e., with a PM larger than 2.3 mas yr"' from the mean 
M 67 motion), even though their colors may have placed them 
near the cluster MS. This particular cut value in the PMs rep- 
resents a compromise between loosing (poorly measured) M 67 
members and including field-type objects. Our final member list 
contains 209 color- and PM-selected objects. In order to mini- 
mize the influence of any uncorrected GD residual, PMs for each 
object were computed with a local sample of members; specifi- 
cally the 25 (at least) closest (r<3'), well-measured cluster stars 
(see Paper I for more details). 

Then we matched our master frame with the sources in the 
Digital Sky Survey (DSS) to compute scale and orientation for 
our master frame. We needed to know scale and orientation with 
an error of ~ 1%. Even if only saturated stars could be matched 



' Note that it was not possible to use the high-probability members 
derived in Paper II as starting reference members, because there are not 
enough stars in common between the two catalogs. Indeed, the PM cat- 
alog of Paper II includes stars in the range magnitude 9<V<21 (with 
reliable PM measurements down to V~19), while our new PM determi- 
nations are in the magnitude range 18<V<26. 
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Fig. 3. Comparison of our derived absolute PMs with values 
from the literature. The three blue points are Hipparcos PM de- 
terminations, while green points are from the Tycho-2 catalog. 
Our M 67 absolute PM is indicated with a red arrow, and the el- 
lipse in its point shows our uncertainty. The other four arrows in 
different colors and types refer to previous values. 

with the DSS catalog, the used sample was good enough for this 
purpose. We divided our displacements for the time baseline be- 
tween the two epochs (10.13 years) to obtain PMs in units of mas 
yr~'. We kept in our final PM list only objects with at least two 
measurements of the displacement, and with PM rms <7 mas 
yr ' in each coordinate. Finally, we corrected our PMs for dif- 
ferential chromatic refraction (DCR) effects, as done in Papers I 
& III, using M 67 white dwarfs and MS stars. We note that DCR 
corrections were always below 0.2 mas yr"' (B - V)"'. 

On the basis of the PM dispersion of members (~0.9 mas 
yr"') with respect to their mean (which reflects our measure- 
ment errors) we estimated the uncertainty on the adopted mem- 
ber reference system to be ±(0.1,0.1) mas yr"'. (Note that our 
estimate of the member dispersion is somehow biased because 
of the 2.3 mas yr"' membership selection criterion we adopted 
in our proper-motion determination. Hence, the dispersion could 
be underestimated. Nevertheless, it still provides a good indi- 
cation of the involved errors.) We were still left with the prob- 
lem of finding the zero point of our PMs - which for now were 
only relative to the cluster's mean motion and not to an absolute 
reference system. Unfortunately, it was not possible to directly 
link our proper motions to the UCAC^J catalog (Zacharias et 
al. |2009t . which includes stars mainly in the 8 to 16 magnitude 
range in a single bandpass between V and R, with no overlap 
with our catalog. 

Background Galaxies can be considered as fixed points in 
the sky, and provide an excellent, and directly-observable abso- 
lute reference system. A visual inspection of the images revealed 
many such galaxies. About 100 of them show point-like nuclei, 
which could be fitted with our PSFs to measure positions and 
PMs (as done in Bedin et al. l2003ll2006l) . 



ihttp://www.usno.navy.mil/USNO/astrometry/optical-IR-prod/ucac/| 



As expected, the errors of galaxy positions are several times 
larger than the typical error of star positions, and depend strongly 
on galaxy morphology. The comparison of the galaxy PMs 
(which should all be the same), suggests that we are underesti- 
mating their errors (which is not the case, for the point sources), 
as a result of a complex combination of seeing and galaxy shape. 
For this reason, we did not use their weighted mean to esti- 
mate the centroid of the galaxy PMs, but adopted an iterative 
cr^ -clipped average instead. 

We started with all the A^'^^lOO galaxies in the sample, for 
which the motion is within 20 mas yr"' from the first-guess mean 
PM of galaxies, and we iteratively estimated their PM dispersion 
cr^ as the 1-D 68.27*-percentile of their distribution around the 
median motion. We then excluded from the sample those galax- 
ies with PM>2.5cr^ from the median galaxies' PM, and we re- 
derived new values of crJJ, median motion and A^*^. We iterated 
this procedure 10 times, noting that the values of median, cr^, 
and A^*-" converged after eight iterations. Final values for cr^ and 
A^*^ are 2.3 mas yr"' and 47, respectively. 

By adopting as the origin of the absolute PM system the fi- 
nal median value of the PM of galaxies, we found for M 67 an 
absolute PM (J2000.0) of 

(jU„cos5,jUi) = (-9.6,-3.7)±(l.l,0.8)masyr"', 

where the uncertainties come from adding in quadrature the un- 
certainties of the centroid of M 67 members to the error of the 
location of the centroid of the galaxies (for each coordinate in- 
dependently). 

The entire sample of ~ 100 galaxies was used to estimate the 
uncertainty associated with the displacement of a single, typi- 
cal galaxy. This was calculated as the 1-D 68.27*-percentile of 
the distribution of the galaxies around their median displacement 
(as computed above). Specifically: (cr'^);j„ cos 5=7.1 mas yr"' and 
{cr'^)^^-5.3 mas yr"'. Because we only used the most consis- 
tent subset of 47 galaxies to derive their relative median mo- 
tion, we associated to the median of their displacements an un- 
certainty that statistically takes into account only those N'^-47 
galaxies. This is done by reducing the error on the single av- 
erage galaxy by the factor 1 / VA'^ - 1. This more-conservative 
approach might result in a overestimate of our internal errors. 
However, considering all the uncertainties involved in the as- 
signment of an error to a galaxy displacement and their limited 
number, this is a preferable approach. 

Figure [T] shows our absolute vector-point diagram (VPD) for 
all the objects in the final list (small dots). M 67 members are 
marked with small crosses. These are the stars within 2.3 mas 
yr"' from the MS stars' mean PM (solid circle). Filled circles 
are visually-confirmed galaxies, and are shown with our esti- 
mated error bars. The best galaxies (selected with the aforemen- 
tioned iterative procedure) are highlighted with red open circles. 
The Figure also shows the 2.5o-'^-5J mas yr"' selection radius 
(dashed circle). The red arrow indicates our estimate of the abso- 
lute PM of M 67. The error bars at the base of the arrow indicate 
the total uncertainty (dominated largely by the estimates in the 
centroid of the relative PM distribution of the background galax- 
ies). In Fig.|2] we used the same symbols as in Fig.[T]to show that 
there are no clear systematic errors in the galaxy PMs as a func- 
tion of the coordinates of the master frame (X, and Y, parallel to 
a and S, respectively). 

The absolute PM of M 67 has been measured by sev- 
eral authors. The first of these determinations comes from van 
Rhijn ( fT922l l. with (//„ cos 5, yu^) = (-7,-h2) + (9,9) masyr"'. 
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Table 1. Non-axisymmetric Galactic model (Pichardo et al. 120031 l2004b . 



Parameter 



Value 



References 



Bar half-length 

Bar axial ratios 

Bar scale lengths 

Bar angle (respect to the Sun) 

Bar mass 

Bar pattern speed (fig) 

Spiral Arms locus 

Spiral Arms pitch angle 

Spiral Arms external limit 

Spiral Arms: exp. with scale-length 

Spiral Arms force contrast 

Spiral Arms pattern speed (Cls ) 



3.1-3.5 kpc 

10:3.8:2.6 

1.7, 0.64, 0.44 kpc 

20° 

10'" Mo 

30-60 km S-' kpc-' 

Bi symmetric (Logthm) 

15.5° 

12 kpc 

2.5 kpc 

~ 10% 

20 km s"' kpc"' 



Gerhard ( l2002t 
FreudetKeich ( I1998t 
Freudenreich ( |1998t 

Gerhard (2002] > 

Debattista et al. ( i2002b 

(•) 

Churchw ell et al . ( l2009t 

Drimmel ( [2000t 

Drimmel ( l2000t 

Disk based 

Patsis et al. ( fTggTT) 

Martos et al. l l2004t 



Note: (•) Weine r & Sellwood ( |1999l >, Fux ( |1999l >, Ibata & Gilmore ( |1995t , Englmaier & 
Gerhard l fT999t . 



Many years later, Murray (11968b computed (jUa cos 6, i^g) - 
(-9.4, -7.0) + (8.0, 2.3) mas yr"'. The more contemporary work 
of Baumgardt et al. (12000il, which made use of the two 
faint Hipparcos stars HIP-43491, HIP-43465 (V-IO), reports: 
QiaC0s6,ns) = (-6 .47, - 6.27) + (1 .29, 1 .01) mas yr"' . Finally, 
Kharchenko et al. (120051) selected 27 M 67 members on the 
basis of their absolute PMs, as derived from the ASCC-2.5 
catalog (Kharchenko 120011 ) — which is based on Hipparcos- 
Tycho family catalogs — and on their location on the V vs. 
B - V color-magnitude diagram. The obtained mean-absolute 
motion of these 27 members is {ij.a cos S, fig) - (-8.31, -4.81) ± 
(0.26, 0.22) mas yr-'. 

Figure |3] shows (in blue) three stars in the field with 
Hipparcos PMs (HIP-43491, HIP-43465 and HIP-43519, with 
error bars), and the stars from the Tycho-2 catalog (in green). 
In the same plot, we mark with a red aiTow our derived M 67 
absolute PM, where the estimated PM error is indicated with an 
elUpse. Note that both the Hipparcos and the Tycho-2 stars are 
far too bright to be measured in our survey. Our M 67 abso- 
lute PM determination is marginally consistent with the bulk of 
Tycho-2 measurements for the objects in the same field. This is 
in line with the expected accuracies for Tycho-2. The same fig- 
ure also shows (dashed arrows) previous determinations of the 
M 67 absolute motion. 

We emphasize that the absolute PM presented here is based 
on direct observations of background galaxies, used to define 
the absolute reference frame. It is a purely differential measure- 
ment, which does not rely, as do previous measurements, on 
a complex registration to the International Celestial Reference 
System through a global network of objects. We end this section 
noticing the good agreement between the absolute PM value of 
M 67 as derived with the bright sources of the Hipparcos cata- 
log (Kharchenko et al. .2005) and that based on the faint, "fix" 
galaxies (this work). 



3. The Galactic orbit of M 67 

With the absolute PM of M 67 given in the previous section, 
its line-of-sight velocity of 33.78±0. 18 kms"' and its heliocen- 
tric distance of 815+50 pc (both from Paper II), we have com- 
puted the Galactic orbit of M 67. We employed a Galactic poten- 



tial that includes the axisymmetric model of Allen & Santillan 
( 11991! ) and the bar and spiral arms models of Pichardo et al. 
(2003, 2004). The axisymmetric background potential of Allen 
& Santillan (.19911) has been scaled to give a rotation velocity 
of 254 km s ' at the solar position, based on the most-recent ra- 
dio astrometry observations by Reid et al. ( 120091 ). We keep the 
original value /?o=8.5 kpc of the solar galactocentric distance 
(Reid et al. ( 12009b give R0-8A+O.6 kpc). The adopted parame- 
ters for the bar and spiral arms and the corresponding references 
are provided in Table [T] The values of the parameters are based 
on recent observations of the Milky Way. The bar model is an 
inhomogeneous ellipsoidal potential that closely approximates 
model S of Freudenreich (fT998l ) from the COBE/DIRBE data 
of the Galactic bar. For the spiral perturbation, Pichardo et al. 
(12003b refined their model until self-consistent orbital solutions 
were found. The orbital self-consistency of the spiral arms was 
tested through the reinforcement of the spiral potential by the 
stellar orbits (Patsis et al. |I99Ib . For an extensive description of 
the models, see Pichardo et al. (l2003l 12004] ). 

Table |2] gives the local standard of rest (LSR) initial ve- 
locity ([/, V, W) and the corresponding cylindrical components 
(n, 0) for three different heliocentric distances of M 67 (cen- 
tral, minimum, and maximum). We used the solar motion 
(U, V, W)=(-10, 5.2, 7.2) km s"' (Dehnen & Binnev [T998] ). U is 
negative toward the Galactic center. In Table [3] we give some 
parameters of the Galactic orbits, corresponding to the three he- 
liocentric distances, computed backwards in time during 1 Gyr 
in the scaled axisymmetric potential. Columns 2 and 3 show the 
minimum and maximum galactocentric distances and Col. 4 the 
maximum z-distance from the Galactic plane; the orbital eccen- 
tricity is given in Col. 5; Cols. 6 to 8 give the azimuthal, radial, 
and vertical periods respectively, and the last two columns the z- 
component of the angular momentum and energy per unit mass. 

Previous computations of the Galactic orbit of M 67 were 
made by Keenan et al. ( 11973b and Allen & Martos ( 11988b in 
other axisymmetric potentials. Keenan et al. (11973b used the 
Galactic models of Schmidt d 1 956b (with an Rq value of 8.2 kpc 
and a circular velocity at the Sun's position yo=216 km s"') and 
of Innanen (1966 ) (with an /?o=10 kpc and a yo=250 km s"'). 
They obtained /?niax=8.8 kpc and 10.7 kpc respectively z,nax= 0.4 
kpc and an eccentricity of 0.1 for both models. Allen & Martos 
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Fig. 4. Meridional orbits using the three heliocentric distances; ii=0.815 kpc (upper panels), 0.765 kpc (middle panels), and 0.865 
kpc (lower panels), computed in the axisymmetric (left panels), barred (central panels), and with bar + spiral arms (right panels) 
scaled Galactic potential. 

Table 2. Local standard of rest initial velocity (U, V, W) and Galactic (11, 0) (kms ') for the three heliocentric distances d (kpc). 



U 



W 



n 







0.815 31.92 ±3.4 
0.765 30.77 + 3.0 
0.865 33.06 ±3.4 



-21.66 ±3.7 -8.71 ±4.3 21.55 ±3.2 233.53 ± 3.8 
-21.03 ±3.4 -6.64 ±3.6 20.98 ± 3.0 234.05 ± 3.5 
-22.28 ±3.9 -10.78 ±4.1 22.12 ±3.4 233.02 ± 4.0 



( |1988t used their galactic model (Allen & Martos |19861 l with values for the solar motion employed. Indeed, the different val- 

7?o=8 kpc and yo=225 km s"'. They obtained /?max=8.7 kpc, ues obtained for /?n,ax result mostly from the scaling of the Sun- 

Zniax=0.48 kpc and an eccentricity of 0.11. In all cases, the or- center distance, as shown by the ratios Ro/Rmax, which are all 

bit of M 67 has a small eccentricity, and the differences found between 1 .07 and 1 .09. The more contemporary, much-improved 

are mostly attributable to the different Galactic parameters and Galactic models and the precise value now available for the ab- 
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Table 3. Orbit parameters in the scaled axisymmetric potential for the three heliocentric distances. 



(kpc) (kpc) (kpc) (kpc) 



P,p P, Pz h E 

(Myr) (Myr) (Myr) (kpc km s"') (lO-km^s^^) 



0.815 


7.65 


9.28 


0.46 


0.096 


207 


146 


76 


2118.3 


-1589.38 


0.765 


7.65 


9.24 


0.42 


0.094 


206 


145 


74 


2114.7 


-1591.50 


0.865 


7.66 


9.33 


0.50 


0.098 


208 


147 


78 


2122.0 


-1587.21 




-5 

X (kpc) 

Fig. 5. Projection of the orbit on the Galactic plane computed in the scaled axisymmetric potential (left), and in the bar + spiral arms 
scaled Galactic potential (right). We adopted an heliocentric distance J=0.815 kpc (full line), 0.765 kpc (dotted line), and 0.865 kpc 
(dashed line). 



solute proper motion of M 67 should result in a more reliable 
orbit for this cluster 

We also computed the Galactic orbit of M 67 in the non- 
axisymmetric potential, first including only the bar, then with 
the bar + spiral arms. The scaled background axisymmetric 
potential was considered. The angular velocity of the bar was 
taken as Qb=60 km s~' kpc"'; the other parameters of the non- 
axisymmetric components are given in Table [T] In the case with 
spiral arms, the mass of these arms was taken as 2.2% the mass 
of the scaled disk component, which is 3% of the mass of the 
original disk. This mass gives a force contrast as listed in Table 

m 

Figure |4] shows the meridional orbits computed in the ax- 
isymmetric potential (left panels), axisymmetric + bar potential 
(central panels) and axisymmetric + bar + spiral arms potential 
(right panels), using the three heliocentric distances. As shown, 
the potential that includes only the Galactic bar gives an orbit 
similar to that obtained with the axisymmetric potential. This is 
because the orbit of M 67 lies far outside the region of the bar 
However, the potential that includes both the bar and spiral arms 
shows a different behavior The orbit is perturbed by the spiral 
arms, mainly in the radial direction. The radial dispersion is not 
very strong in M 67, but distorts what would be a box orbit. Thus, 
a moderate spiral potential has important effects in the kinemat- 



ics of orbits near the Galactic plane, as is the case in M 67, and 
in general, as is the case for the solar neighborhood stars (Antoja 
et al. l2009l l. This result holds for the allowed variations in proper 
motion, radial velocity and distance, because the most important 
parameter that affects the orbit is the mass of the spiral arms. 

In left panel of Fig. |5]we show the projection of the orbit on 
the Galactic plane, computed in the axisymmetric potential. At 
the scale shown in this figure there is no appreciable difference 
between the orbits using the three heliocentric distances. The 
right panel shows the corresponding orbits in the axisymmetric 
+ bar + spiral arms potential. There is a slight difference in the 
azimuthal behavior in both figures. 

We have also computed the tidal radius of M 67 in the ax- 
isymmetric + bar + spiral arms scaled potential. With a total 
mass in M 67 of 279 Mp, lis ted by Piskunov et al. (520081), using 
King's equation (King [T962] l we obtained a mean tidal radius of 
7. 1 pc. With the alternative equation (1) in Allen et al. (12006b the 
tidal radius is 8 pc. Both results are near the 9.6 pc value listed 
by Piskunov et al. (2008). 
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